function [ trendmas, good, baddays ] = TrendDef( times, mas )
%TrendDef Trend defenition + good data marker
        trendmas=csaps(times,mas,1e-12,times);
        good=abs(mas-trendmas)<4*sqrt(trendmas);
        good=good==1;
        
        trendmas=csaps(times(good==1),mas(good==1),1e-12,times);
        good=abs(mas-trendmas)<4*sqrt(trendmas);
        good=good==1;
%         pol=polyfit(times,mas,5);
%         trendmas=polyval(pol,times);
%         good=abs(mas-trendmas)<4*sqrt(trendmas);
%         good=good==1;
        baddays=sum(good)./length(good);
%         clf
%         plot(times(good==1),mas(good==1),'-');
%         hold on
%         plot(times,trendmas,'r-');
%         pause(0.1)
end

